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■j^' ■ Abstract 

' This paper characterizes the effect of finite rate channel state feedback on the sum rate of a muhi-access 

. multiple-input multiple-output (MIMO) system. We propose to control the users jointly, specifically, we first choose 



the users jointly and then select the corresponding beamforming vectors jointly. To quantify the sum rate, this paper 
'. introduces the composite Grassmann manifold and the composite Grassmann matrix. By characterizing the distortion 

rate function on the composite Grassmann manifold and calculating the logdet function of a random composite 
Grassmann matrix, a good sum rate approximation is derived. According to the distortion rate function on the 
' composite Grassmann manifold, the loss due to finite beamforming decreases exponentially as the feedback bits 

on beamforming increases. 
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I. Introduction 

^ ' This paper considers the uplink of a cellular system with one base station and multiple users, where 
both the base station and each user are equipped with multiple antennas. Multiple antenna systems, also 
C\l • known as multiple-input multiple-output (MIMO) systems, provide significant benefit over single antenna 
^ ' systems in terms of either higher spectral efficiency or better reliability. For the uplink of a cellular system, 
it is reasonable to assume that the base station has the full knowledge about the uplink channel while the 
^ users has partial information about the uplink channel through a feedback link from the base station. In 
O practice, it is also reasonable to assume that the feedback link is rate limited. 

j> ! The purpose of this paper is to quantify the effect of the finite rate channel state feedback on the sum 
rate. The effect of finite rate feedback on single user MIMO systems are well studied. MIMO systems 



with only one on-beam are considered in [1] and [2] while systems with multiple on-beams are discussed 



in [3]-[8]. In the recent works [7] and [8], the effect of finite rate feedback is accurately quantified 
by characterizing the distortion rate function in the Grassmann manifold. For multi-access systems, the 
throughput capacity region is characterized in [9] with the assumption that each user has only one antenna 
and the full channel information is available to all users. 

To characterize the feedback gain, we propose to control the users jointly. An simple extension of [10] 
can show that the optimal strategy is to select the covariance matrices of the transmit signals of the users 
jointly. It is different from the current systems where the base station controls the users individually. 
The gain of joint control over individual control is analogous to that of vector quantization over scalar 
quantization. However, it is difficult to either implement or analyze the fully joint control. For simplicity, 
this paper proposes a suboptimal strategy employing power on/off strategy, where we first choose the on- 
users jointly and then select the beamforming vectors jointly. The effect of user choice can be analyzed 
by extreme order statistics. To quantify the effect of beamforming, the composite Grassmann manifold 
is introduced in this paper. By characterizing the distortion rate function on the composite Grassmann 
manifold and calculating the logdet function of a random composite Grassmann matrix, a good sum rate 
approximation is derived. According to the distortion rate function on the composite Grassmann manifold, 
the loss of finite beamforming decreases exponentially as the feedback bits on beamforming increases. 
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II. System Model 

Assume that there are Lr antennas at the base station and N users communicating with the base 
station. Assume that the user i has » antennas 1 < i < A^. In this paper, we let Lr,i = Ltj — Lt for 
^ ^ i, j ^ N. The signal transmission model is 

TV 



Y = ^HiTi + W, 



where Y e C^^^^ is the received signal at the base station, Hj e c^axLt channel state matrix for 
user i, Tj is the transmitted Gaussian signal vector for user i and W G C^^^^ is the additive Gaussian 
noise vector with zero mean and covariance matrix 1^^. In this paper, we assume the Rayleigh fading 
channel model, i.e., the entries of Hj are independent and identically distributed (i.i.d.) circularly symmetric 
complex Gaussian variables with zero mean and unit variance (CJ^ {0, 1)) and Hj's are independent for 
each channel use. 

We assume that there exists a common feedback link from the base station to all the users. At the 
beginning of each channel use, the channel states H/s are perfectly estimated at the receiver. A message, 
which is a function of the channel state, is sent back to all users through a feedback channel. The 
feedback is error-free and rate limited. The feedback directs the users to choose their Gaussian signal 
covariance matrices. In multi-access system, users are uncoordinated. It is reasonable to assume that 

t t 

T[ • • • T]y be the overall transmitted Gaussian signal for all users and 



E 



T,:T!. 



= 0. Let T 



S = E [TT^] be the overall signal covariance matrix. Then S is an NLt x NLt block diagonal matrix 



TT' 



Assume there is a covariance 



whose i^^ diagonal block is the Lt x Lt covariance matrix E 
matrix codebook B-z = {Si, • • • , "^Kb} declared to both the base station and the users, where each 
is a proper overall signal covariance matrix and fCg is the size of the codebook. Let H = [H1H2 ■ ■ ■ Hat] 
be the overall channel state matrix. The feedback function 99 is a mapping from {He j^Lrxnlt^ ^^^q 
the index set {1, • • • , Kb}. Subjected to the finite rate feedback constraint 

l^sl < Kb 

and the average transmission power constraint 

Eh [tr(S^(H))] <p, 
we are interested in characterizing the sum rate 

max max Er [log IIl^ + HS^(H)H"f|l . (1) 

Since the variance of the Gaussian noise is normalized, the average power constraint p is also the average 
received signal-to-noise ratio (SNR). 

III. Mathematical Preliminary Results 

For compositional clarity, this section assembles the useful mathematical results that we derive for later 
analysis. Due to the space limit, we omit all the proofs. 

A. Extreme Order Statistics for Chi-Square Random Variable 

Let Xj = \ where hij 1 < j < L, 1 < i < n are i.i.d. circularly symmetric complex 

Gaussian variables with zero mean and unit variance. Let us rearrange these i.i.d. chi-square random 
variables Xi, ■ ■ ■ , X„ into a nondecreasing sequence Xj^ < Xi^ < • • • < Xj„. Let n approach infinity, the 

Yl''k=i ^in-k+i where I is a fixed positive integer. 



following theorem gives a formula for E 
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Theorem 1: Let X — Y^^=i l^jl^ where hj ~ CAf {0, 1). Denote the distribution function of X by 
Fx (x). Then for any fixed positive integer I, 



lim E 

n— *+oo 



where a„ is the solution of 



and 



/it 



k=l 



mi\x: 1-Fx{x)<-}, 

EL-l L-i i 
i=0 i\ "n 



xde~ 



0.577216- •• . 



Although this theorem is for asymptotically large n, it gives an accurate approximation when < / -C n. 



B. Conditioned Eigenvalues of the Wishart Matrix 

Let H e L"^"* be a random nxm matrix whose entries are i.i.d. Gaussian random variables with zero 
mean and unit variance, where L is either M or C and m < n w.l.o.g.. The random matrix W = H^H is 
Wishart distributed and its distribution is denoted by Wm (n^lm). 

For a W ~ {n, I^), the following proposition shows that conditioned on the trace, the conditional 
expectation of a specific eigenvalue of W is proportional to the condition with a ratio independent of that 
condition. 

Proposition 1: Let W ~ Wm (n, I„i) where n> m. List the ordered eigenvalues of W as Ai > A2 > 
• • • > Am > 0. Then conditioned on the trace of W, i.e., ^™ ^ \ = c where c > 0, the ratio between the 
conditional expectation of Aj and the condition c is a constant Q independent of c, i.e.. 



E 



Ail Aj 



1=1 



where 



/^A,=l ^illJLlA/ 



(n-m+l)-l 



A„.(A)i^nr=i^^.- 



Al>-->A„ 



|(n-m+l)-l 



^1 ^ " rifATTj, 



/5 = 1 if L = R or /3 = 2 if L = C, and |Am (A)| = {\ - Xj). 

In general, it is not easy to calculate the constant Q I < i < m. Fortunately, the constants can be well 
approximated by asymptotics. Due to the space limit, we only present the asymptotic formula for Ci in 
the following proposition. 

Proposition 2: Let the random matrix W ~ (n, 1^) where n > m. Define y = ^. Then the 
asymptotic approximation gives 



E 



1 

71 



TT — a + - sin (2a) 



where a satisfies 



and 



1 

m 



TT — a ^ sin (a) + ^-^6'„ 



TT — a — sin (a)] 



if 2/ < 1 

if y = i 



9y — tan ^ 



y/ysin (g) 
- ^ cos (a) J ■ 
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C. The Grassmann Manifold and the Composite Grassmann Manifold 

The Grassmann manifold is the geometric object relevant to the beamforming quantization analysis. 
The Grassmann manifold '3n,m (L) is the set of m-dimensional planes (passing through the origin) in 
Euclidean n-space L". A generator matrix P G L"^*" for an m-plane P E Qn,m (L) is the matrix whose 
columns are orthonormal and span P. The generator matrix is not unique. That is, if P generates P then 
PU also generates P for any mxm orthogonal/unitary matrix U (w.r.t. L = M/C respectively) [11]. The 
chordal distance between two m-planes Pi, P2 E Gn,m (L) can be defined by their generator matrices Pi 



PiPj - P2PI 



[11]. The uniform distribution on Gn,m(^) with density 



1 

satisfies fp (Pi) = fp (P2) for arbitrary Pi, P2 e gn,m (L) [12]. 



and P2 via dc{Pi,P2 
function fp 

For quantizations on ^„ ^ (L), the corresponding distortion rate function has been characterized [7]. A 
quantization q on ^„ m (L) is a mapping from Qn,p (L) to a subset of Gn,p (L), which is typically called 
a code C, i.e., q : Qn,p (L) C. Define the distortion metric as the squared chordal distance. Then the 
distortion associated with a quantization q is 



D^Ec 



[dliQ,qm] 



where the source Q is randomly distributed in Qn,m (L)- Assume that the source Q is uniformly distributed 
in Qn,p (L) . For any given code C, the optimal quantization to minimize the distortion i^ 

q (Q) = arg mm 4 {P, Q) ■ 



The distortion associated with this quantization is 

D{C)=Eq 



min dl (P, Q) 



For a given code size K where K is a positive integer, the distortion rate function i^ 

D* (K) = inf D (C) . 

C:\C\=K 

In [8], we derive a lower bound and an upper bound for L = C 

t 1 logo if r (t) 1 logo if 



t+1 



where t = m(n — m), 



— TT™ {n—iy. 
f! lli=l (m-iy. 
j_ y-rn-rr 
t\ lli=l 



t 



if 1 < m < f 



, ^n-m (n-iy if n < ^ < ^ ' 
(n—m—t)\ 2 — — 



and the symbol < denotes the main order inequality, f (K) < g (K) if lim 4^ < 1. 

To treat multi-access MIMO systems, we define the composite Grassmann manifold. The /c-composite 
Grassmann manifold Qnjn (L) is a Cartesian product of k Qn,m (IL)'s. Denote P^^^ an element in Qnjn (L). 

p(^) = (Pi,--- ,P,) 



where Pi G Gn,m (L) I < i < k. For any Pi'''', P2'''' G Qn^L (L), we define the chordal distance between 
them 



d^{p['\p^'^)= j^dl{P,,,P,, 
\ i=i 



'The ties, i.e., the case that 3Pi, P2 G C such that dc {Pi, Q) = min dc {P, Q) = dc {P2, Q), are broken arbitrarily because the probability 
of ties is zero. 

^The standard definition of the distortion rate function is a function of the code rate defined by log2 K. The definition in this paper is 
equivalent to the standard one. 
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where P^'^^ — (Pi,i, • • • , Pi,k) and P2*^^ — (-^2,1, • • • , P2,k)- It is easy to verify that the chordal distance 
on Q^li (L) is well defined. 

This paper characterizes the distortion rate function for quantizations on Qn,m (L). Define the distortion 
metric on Qnjn (L) as the square chordal distance on it. Assume a uniformly distributed source Q^''^ in 
Qn,L (L). The following theorem characterizes the distortion rate function for quantizations on Qn,L (C). 

Theorem 2: The distortion rate function on Qn li (C) is upper bounded and lower bounded by 



1 _ / 1 N . J_ 

kt 



A;t + 1 Vr(H + 1) ' y ~ ^ ' ^ kt \T{kt + l) 

where t = m{n — m), 

ifl < m < f 

=1 (n— m— ij! 2 — — 

and the symbol < denotes the main order inequality, f (K) < g (K) if lim 4^ < 1. 

It is noteworthy that the upper bound is derived by computing the average distortion over the ensemble 
of random codes. In practice, we often use the upper bound as an approximation to the actual distortion 
rate function. 




D. Composite Grassmann Matrix 

Roughly speaking, a composite Grassmann matrix is the generator matrix for an element in Qnjn (L). 
Let P^'^) = (Pi, ■ ■ ■ ,Pfc) G giSn (L). The composite matrix p('=) generating p(*^) is P^'^) = [Pi ■ ■ ■ P^] 
where Pi, ■ ■ ■ , P^ are the generator matrices for Pi, ■ ■ ■ , Pfc respectively. Since the generator matrix for 
a plane in the Grassmann manifold is not unique, the composite Grassmann matrix generating P^'^) is not 
unique either. Let P^*^^ be a generator matrix for P^'^\ The matrix p('^)U^'^), where U^'^^ is the arbitrary 
km X km block diagonal matrix whose k diagonal blocks are m x m orthogonal/unitary matrices (w.r.t. 
L = M/C respectively), also generates P^'^K In this paper, the set of composite Grassmann matrices for 



gL% (L) is denoted by M\^^,L (L). 

For a random composite Grassmann matrix P^'^), the following theorem bounds E [log |l + cP^'^^^P*^'^) |] . 
Theorem 3: Let P^'^^ e -^ni (^) be uniformly distributed. For any positive constant c. 



(k) 



E 



H 



log 



Ife+-HtH 

n 



< 



Epk [logllfc + cP^^'^tpe^^l] 
< logEp.[|l, + cP('=)tp«|], 
where H G L"^^ has i.i.d. Gaussian entries with zero mean and unit variance. 

In the above theorem, both bounds can be computed explicitly. In [13], we derive an asymptotic formula 
to approximate the lower bound. Let n and k approach infinity simultaneously with fixed ratio. 



lim 



(n,fe)-»+oomin (n, k) 



E 



H 



log 



Ik + -H^H 



n 



, , , , , u (1 —v) log (1 — ur) 

= log {w) - log a ^ ^^-^ 

r y 



1 A min(?i,A;) A 

where y — — 7-^, r = 



mm[n,k)-c 



, W 



A 1 



l + y + a+ \/{l + y + a) — % and u = 



^ ^l+y + a — \J {l-\-y -\- of — . Formulas for the upper bound are also derived in this paper. Due 

to the space limit, we only present the formulas for 1 < A; < 5. The expectation Ep/t [|lfe + cP*^'^)^P^'^^|] 
can be calculated by 
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k=l 1 + c; 
k=2 (1 + c? - c^i- 



n ' 



k=3 (l + cf -c2(l+c)! + c3A; 



k=4 (l + c)^-c2(l + cf f + cMl + c)^-c^(;^-^); and 

k=5 (l + cf-cMl + cf i^ + c3(l + cf f -cni + c)(f -i)+c^(i-f). 

IV. The Suboptimal Strategy and the Sum Rate 

This section is devoted to calculate the sum rate of a multi-access MIMO system with finite rate 
feedback. The computation of the sum rate ([U) involves two correlated optimization problems: one is 
with respect to the feedback function Lp and the other optimization is over all possible covariance matrix 
codebooks. The direct calculation of ([B is difficult. 

To reduce the complexity, we propose a suboptimal strategy to control the users jointly. Specifically, 
we first choose the on-users jointly and then select the corresponding beamforming vectors jointly. It is 
different from the current system where users are controlled individually. 

The assumptions for transmission are as follows. 

Tl) Power on/off strategy. In power on/off strategy, the user i's covariance matrix is of the form = 
-PouQj QL where Pon is a fixed positive constant to denote on-power and Q,; is the beamforming 
matrix for user i. Denote each column of Qj an on-beam and the number of the columns of Qj 
by li, then Q,|Qi = I;, where < k < Lt and k = is for the case that the user i is off. This 
assumption is motivated by the fact that power on/off strategy is near-optimal for single user 
MIMO systems [7]. 

T2) At most one on-beam per user. This assumption implies either /j = or /j = 1. It is proposed 

so that each user has larger probability to be turned on. 
T3) Constant number of on-beams for a given SNR. Let / = J2iLi ^« be the total number of on-beams. 

we assume that I is a constant independent of the specific channel realization for a given SNR. 

This assumption is motivated by the fact that constant number of on-beams is near optimal for 

single user systems [7]. It will be validated for multi-access systems in later analysis. 
The feedback is described as below. 

Fl) User selection criterion. Assume that / users will be turned on. We choose the / users with the 
largest channel state Frobenius norms, i.e. ||Hj^ || > ||Hj|| for all i ^ {ij : I < j < 1} where ||-|| 
is the Frobenius norm and ii, - ■ ■ ,ii are the users chosen to be on (on-users). 
According to this user selection criterion, the feedback contains two parts, one of which indicates the 
/ on-users and the other of which is for beamforming. Let ii,--- ,ii be the on-users and bi,--- , b; 
be the beamforming vectors for those users. Then B = [bi ■ ■ ■ b;] G A^^^ i (C) where M-^^i (C) 
is the set of composite Grassmann matrix (Section IIII-D|) . Denote the beamforming codebook B = 
|Bfc : Bjt G M-Lt 1 (^) ) 1 ^ ^ ^ l-^lj- Then the overall feedback codebook is the Cartesian product of 
the set of on-users {(ii, ■ • ■ and the beamforming codebook B. Let Hj^ be the channel state matrix 
for the user ij and i be the right singular vector corresponding to the largest singular value of Hj^ . 
Define V = [vi j • ■ ■ i]. The beamforming feedback function is defined as the following. 

F2) Beamforming feedback function. 



99([H,,---H,J) ^ arg min (V, B^) 

l<fc<|B| 

I 



arg max vj-i^b^ ,,- 



l<k<\B\ .^^ 



(2) 



where (V, B^) denotes the chordal distance between the elements in the composite Grassmann 
manifold ^ (C) generated by V and Bfc, and hkj is the j*^ column of the k^^ beamforming 
matrix Bi, G -B. 



7 



The feedback assumptions Fl and F2 will be validated in the later analysis. 

The above assumptions define a suboptimal strategy for multi-access MIMO systems. The key point 
is that the user choice is independent of the channel directions and the beamforming is independent of 
the channel strengths (norms). In this way, the effect of user choice and beamforming can be studied 
separately. Before diving into the general analysis, we discuss a special case, antenna selection, to get 
some intuition. 



A. Antenna Selection 

The system model for antenna selection is 



Y=J2 h.T. + W, 



i=l 



where hj is the i^^ column of the overall channel state matrix H. For each channel realization H, we 
simply choose / antennas zi, ■ ■ ■ ,ii such that ||hjj| > ||hj|| for all i ^ {ij : 1 < j < /}. Here, we actually 
do not require one on-beam per on-user (Assumption T2). Write h^. = nj^j where rij is the Frobenius 
norm of hj^ and is the unit vector to present the direction of hj^ . Define S = [^i • ■ ■ ^/]. We have the 
following upper bound on the sum rate. 



E- 



H 



log 



= Eh 
< Es 

= E= 



log 
log 

log 



diag [nf, 



n; 



1/ + ^E„2 [diag [nl, 



1^ + 



(3) 



where the inequality follows from the concavity of log|-| function and the fact that ra^'s and H are 



1 1 2 

independent. Noting that ||hj|| 's are i.i.d. chi-square random variables, an accurate approximation to 



En2 Yl'j=i tie obtained for / ^ by applying the asymptotic extreme order statistics in Theorem 

[T] On the other hand, it can be proved that ^/s are independent and uniformly distributed unit vectors. 
Regarding S as a Cartesian product of ^j's, 3 is also uniformly distributed in A^^^^ i, the set of composite 
Grassmann matrix. According to the results in Section UlI-DI for E [log |l + cS^^S|], the upper bound of 
the sum rate dH) can be characterized. Simulation show that the upper bound Q is tight. The sum rate of 
antenna selection is then approximately characterized. 



B. General Beamforming 

With the assumptions Tl-3, Fl and F2, the signal model for the general beamforming is 

I 

Y = 5^H,^(H),T, + W, 

i=i 

where b(^(H),j is the column of the feedback beamforming matrix B^(h) € B. For notational conve- 
nience, we denote b(p(H),j by b* and the equivalent channel vector Hj^.b* by hj. Let rij be the Frobenius 



8 



norm of hj, be the unit vector presenting the direction of hj and S = [^i ■ ■ ■ ^/]. Then the sum rate is 
given by 



E 



H 



E 



log 



H 



log 



1; + ydiag [ni, ■ ■ ■ ,nf] H'^H 



It can be proved that ^/s are uniformly distributed and independent of n/s. Denote the singular 
value decomposition of Hj^ by UjAjVj. After beamforming, the equivalent channel vector for user 

ij is hj = Uj ^AjVjb*j = Uj^jUj, where is the direction of the vector AjVjb*. Since the user 
choice is only dependent on Aj's and the beamforming matrix selection is only relevant to V^'s, U/s are 
independent and uniformly distributed. According to [14, Thm. 6.1], = Vj^j is uniformly distributed 
and independent of rij's. Thus, similar to (|3]), the sum rate of general beamforming can be upper bounded 
by 



Es log 



EL. 



(4) 



where H = [^i ■ ■ - ^j]. 

It is more involved to calculate E 
rt 



Ej=i ■ Let Xj^k I < k < Lt he the ordered eigenvalues of 
H.],H.i^ such that Aj^i > Xj^2 > ■ ■ ■ > Xj^l^ > 0. Let ^ be the right singular vector of Hj^ corresponding 
to the k^^ largest singular value a/A^. Then 



E [r^] 



E|||H,^b*| 



E 



b*^H^ H b* 



E 



.k=l 

k=l 



(5) 



where the last equality follows from the fact that the beamforming is independent of the channel norms, 
i.e., J2k=i ^j,k^- The E [Aj J's can be calculated by 



E[A, 



E 



E 



A 



CfcE 



H,- 



(6) 



where the last equality is a direct application of Proposition [T] in Section Ull-Bi To evaluate E 
we need the following proposition. 



<kbj 



Proposition 3: Consider the beamforming feedback function in (|2l). Define 7 = E 



Then E 



and E 



(1 - ^) / {Lt - 1) for all 1 < j < / and 2 < A; < Ly. 



Apply this proposition and substitute Q into ([5]). After some elementary manipulations, we have 



E 



■ / 
.i=i 



Ci7 (l-Ci)(/-7) 



T 



(7) 



H,; 



and (i 



Theorem [T] and Proposition [2] provide asymptotic formulas to approximate Yl]=i E 
respectively. Define K = \B\ the size of the beamforming codebook. The maximum 7 achievable 7sup 
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is a function of K. According to the distortion rate function D* (K) for quantizations on the composite 
Grassmann manifold Qf^^ ^ (C), 

7sup = sup = 1 - D* {K) . 



B: \B\<K 



Substitute 7sup into ([7]). The expectation E Yl]=i 



can be calculated as a function of K. 



Finally, substituting the value of E Yl]=i 



into dH) and employing the results in Section IIII-DI for 

E [log 1 1 + cH1^H|] , the upper bound of the sum rate (H)) can be characterized. Simulations show that this 
upper bound is tight. The sum rate is therefore approximately characterized. 

C. The Effect of Finite Rate Feedback 

The above analysis characterizes the effect of finite rate channel state feedback. The upper bound (H]) 



shows that the effect of feedback is quantified by E 
choice and beamforming can be analyzed separately" 

According to dV]), the effects of user choice is reflected by Yl\=i ^ 



. Formula (|5]) shows that the effect of user 



sum rate requires to maximize Yl]=i E 
is validated. Furthermore, the term Xl'ti ^ 



H,; 



Maximization of the 



and thus the user selection criterion (Assumption Fl) 



is an increasing function of the number of users 
(Refer to Section UlI-AI) . The more users the system"" has, the larger the sum rate is. 

The effect of beamforming can be analyzed according to dV]). Define K = \B\ and i?fb = \0g2K. 

Assume that i?fb is large so that (/ — 7) ^ 7. Then approximately, E Yl]=i ^"j proportional to 7. 
The beamforming feedback function should maximize 7 and Assumption F2 is therefore verified. Denote 
I ~ 7sup = D* (K) the beamforming loss. From the distortion rate function on the ^ (C), / — 7sup is a 
exponentially decreasing function of Rfh/l (Lt — !)• We expect that a few feedback bits on beamforming 
could have large gain while more feedback bits wouldn't gain much further. 

The assumption T3 about the constant number of on-beams can be validated as well. Assume that both 
the number of users and the feedback bits on beamforming i?fb are large. Because of the user choice 
and beamforming, the quantities 1 < j < I are relatively "stable", i.e., the fluctuations of n^'s are 
relatively small. It is reasonable to assume constant number of on-beams for multi-access system. 

The antenna selection can be viewed as a special case of general beamforming where the beamforming 



vector is always a column of the identity matrix. For general beamforming, log2 



+ i?fb feedback 



/ NLr \ ( N \ 

bits are needed. For antenna selection, there are logg I ^ 1 ~ log2 I ^ j ^"^^2 feedback bits 

needed. Since antenna selection does not assume one on-beam per on-user (Assumption T2), it is expected 
that the sum rate of antenna selection is close to but better than that of general beamforming with the 
same / and -Rfb = / log2 Lt- The improvement is due to the extra freedom the antenna selection has. 

D. Simulation 

The sum rates of antenna selection and general beamforming are given in Fig. [Hand Fig. |2] respectively. 
Simulations show that the upper bound (H]) (solid lines) is tight. Note that the upper bound dH) is of the 
form E [log |l + cS^S|] . Theoretical analysis (Theorem [3]) gives an upper bound (plus markers) and a 
lower bound ('x' markers) on dll). Simulations show that these theoretical approximations are accurate. 

Fig I2] also depicts the gain of beamforming. The sum rate by finite rate beamforming feedback (circles) 
is compared to that of perfect beamforming (dash-dot lines). Simulation shows that with several feedback 
bits on beamforming, the corresponding sum rate is close to that of perfect beamforming. As a special case 
of general beamforming, antenna selection is shown to be similar to but better than general beamforming 
with the same / and i?fb = / log2 Lt. 
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V. Conclusion 

This paper proposes a strategy where users are controlled jointly. The effect of user choice is analyzed 
by extreme order statistics and the effect of beamforming is quantified by the distortion rate function 
in the composite Grassmann manifold. By characterizing the distortion rate function on the composite 
Grassmann manifold and calculating the logdet function of a random composite Grassmann matrix, a 
good sum rate approximation is derived. 
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